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ABSTRACT 

We present a spatially resolved spectral analysis of full ASCA observations of the remnant of the 
supernova of 1006 AD. This remnant shows both nonthermal X-ray emission from bright limbs, gen- 
erally interpreted as synchrotron emission from the loss-steepened tail of the nonthermal electron 
population also responsible for radio emission, and thermal emission from elsewhere in the remnant. 
In earlier work, we showed that the spatially integrated spectrum was well described by a theoretical 
synchrotron model in which shock acceleration of electrons was limited by escape, in combination 
with thermal models indicating high levels of iron from ejecta. Here we use new spatially resolved 
subsets of the earlier theoretical nonthermal models for the analysis. We find that emission from 
the bright limbs remains well described by those models, and refine the values for the characteristic 
break frequency. We show that differences between the northeast and southwest nonthermal limbs 
are small, too small to account easily for the presence of the northeast limb, but not the southwest, 
in TeV gamma-rays. Comparison of spectra of the nonthermal limbs and other regions confirms that 
simple cylindrically symmetric nonthermal models cannot describe the emission, and we put limits on 
nonthermal contributions to emission from the center and the northwest and southeast limbs. We can 
rule out solar- abundance models in all regions, finding evidence for elevated abundances. However, 
more sophisticated models will be required to accurately characterize these abundances. 
Subject headings: 



I. INTRODUCTION 

A rapidly growing number of galactic remnants ex- 
hibit nonthermal (and non-pulsar-related) emission in 
the X-ray band. Nonthermal emission can dom- 
inate over much or all of the remna nt, as is 
the case with SN 1006 fDver et al." "2000), G266.2- 
1.2 [RX J08 52.0-46221 (Slane et al. 2001), G347.3-0.5 
llSlane et al.l Qj399; Uchivama. Aharonian. & Takahashi 
I2003() . and AX J1843.8-0352, (,Ueno et aL.2n0 3'l. Alter- 
natively an SNR can have signatures of nonthermal emis- 
sion along with X-ray spectral lines, for example RCW 86 
(flho et al. 20 03ii'), Cassiope ia A (Vink & Lamine 2003), 
and Tycho, ( Hwang et alj[2002|) '. Finally some SNRs 
(Cassiopeia A, Kepler, Tycho, SN 1006 and RCW 86) 
reveal nonthermal emission in the hard- X-ray sensitivity 
of the RXTE PGA (up to 60 keV), suggesting nonther- 
mal emission may be common among most young SNRs 
HAllen. Gotthelf PetrellTofl^. 

In most of these instances, a strong case can be made 
that the nonthermal emission is synchrotron radiation. 
The total number of SNRs with suspected nonthermal 
emission is now in the double digits and we expect, as 
new X-ray instruments improve in spectral and spatial 
resolution, many more SNRs will be found with vary- 
ing amounts of X-ray synchrotron emission, previously 
undetected among thermal continuum and line emission. 
In the past this emission was missed because few anal- 
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yses included models for nonthermal emission. Even of 
those that do, many do not take advantage of this unique 
opportunity to extract information about particle accel- 
eration from the emission of ultra-relativistic electrons. 
We know that the spectrum extrapolated from radio fre- 
quencies must roll off before X-ray energies to avoid ex- 
ceeding X-ray flux measurements (Reynolds & Keohane 
1999; Hendrick & Reynolds 2001). Therefore observa- 
tions of X-ray synchrotron emission lie in the regime 
where the particle spectrum is changing rapidly. Shock- 
acceleration models that can reproduce this drop-off may 
be able to provide information about the remnant age, 
radiative losses of electrons or the spectrum of magneto- 
hydrodynamic (MHD) waves near the shock. 

In Dyer et al. (2000; hereafter Paper I) we demon- 
strated that the total-ffux spectrum from SN 1006 was 
well fit by a combination of new synchrotron models and 
plane-shock thermal models, and that the fit was an im- 
provement over previous fits both due to the robustness 
of the synchrotron models (which use radio measure- 
ments to constrain the nonthermal emission) and due 
to the new accuracy of the thermal models, which, with 
the assistance of the synchrotron model, detected for the 
first time half a solar mass of iron. However, these mod- 
els were fit to spectra summed over the entire remnant 
and it is well known that the remnant spectra can vary 
significantly across t he face of a remnant. In particular 
iKovama et al.l l)1995tl settled a controversy about the X- 
ray spectrum of SN 1006 by exhibiting a difference in 
spectra between the limb and center. 
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One of the most compelling reasons for spatially re- 
solved studies of SN 1006 is the nonaxisymmetric TeV 7- 
ray emission from SN 1006 l)Tanimori et al.ll998|) . In Pa- 
per I we interpreted this emission as inverse- Compton up- 
scattering of cosmic-microwave-background photons and 
derived a preshock magnetic field of ~3 /^G, implying a 
mean field in the remnant of about 9 fiG. Since then, full 
analysis of the 7-ray spectrum yielded a magnetic field 
of 4 /iG ijTanimori et aljr2001f) . in reasonable agreement 
with our previous results. Inverse-Compton is not the 
only possible emission mechanism to explain the 7- rays 
- recent work l)Berezhko. Ksenofontov. & Volkll2002t) at- 
tributes the TeV emission from SN 1006 to decay of 7r° 
mesons produced by inelastic collisions between cosmic- 
ray protons and thermal gas. However, neither interpre- 
tation explains why 7-rays are seen only in the northeast 
(NE) - only spatially resolved studies could hope to an- 
swer this question. 

The clear symmetry of SN 1006 about a northwest- 
southeast axis suggests that generalizations of the cylin- 
drically symmetric model used in Paper I (SRESC) is a 
reasonable starting point for such studies. While the 
symmetric synchrotron model is an oversimplification, 
the results from the whole remnant fits (Paper I) were 
very encouraging, suggesting that we extend the model 
to find its limits. We hope to answer a set of related 
questions: Is the spatial distribution of emission from 
SN 1006 consistent with an axially symmetric model? 
How exact is the bilateral symmetry, and can any hint 
be found in the X-rays for the gross asymmetry shown by 
the TeV 7-ray emission? In the regions where nonther- 
mal emission dominates can thermal emission be charac- 
terized by simple models? 

There has bee n much recent work on SN 1006. 
iVink et aTl l)200 0fl analyzed integrated spectra from Bep- 
poSAX of SN 1006 with a new version of SPEX designed 
to simultaneously fit regions which overlap in the point 
spread function. They found an adequate two temper- 
ature fit to the thermal emission and supersolar abun- 
dances. New high resolution Chandra observations of 
the NW and NE provide the best spatial resolution im- 
ages to date and several groups have been exploring the 
fine details. iLong et "aTl 1)200,*^ analyzed the first two 
Chandra poi ntings and foun d no evidence for the halo 
predicted bv iRevnoldsl l)1996j) . The radio and X-ray fea- 
tures were found to be perfectly correlated in the NE al- 
though above 0.8 keV X-ray limb brightening is more pro- 
nounced. Clumps of thermal X-ray material were found 
interior to the shock in t he NW and NE with super-solar 
abundances in the NW. iBamba et alJ |)200.1D carried out 
spatial-spectral fits of cross sections of the Chandra ob- 
servation of the sharp fine filaments in the NE. Using 
an energy cut to separate thermal and nonthermal emis- 
sion they found the nonthermal emission to be narrower 
than the thermal emission. From thermal fits they found 
super-solar abundances, including iron, and a profile con- 
sistent with Sedov dynamics. More recent Chandra ob- 
servations by Hughes use short exposures to image the 
entire SNR. 

Here we attempt to extract the maximum useful in- 
formation from all ASCA observations of SN 1006. We 
present spatial subsets of the SRESC synchrotron model 
of Paper I, and apply them to describe the nonther- 
mal emission. We believe that this study represents the 



fullest use of the substantial amount of observing time 
obtained with the ASCA satellite on SN 1006, and that 
further advances will require both better data and better 
thermal and nonthermal models. 

2. THE CONTEXT OF THIS WORK 

Much ground has been covered in the process of moving 
from phenomenological power laws to the development of 
a synchrotron model appropriate for regions of SN 1006. 
[Revnolds & Keohanc ( 1999) used a maximally curved 
model, SRCUT to find upper limits for synchrotron emis- 
sion in Galactic remnants. The model was fit, ignoring 
evidence for thermal emission, as if all X-ray emission 
were synchrotron. Since the model had maximal curva- 
ture, this procedure produced the highest electron en- 
ergy at which the lower-energy power-law could begin 
to steepen and still not exceed observed X-rays, placing 
a solid upper limit on the energy to which SNRs could 
accelerate electrons with the same slope as at radio emit- 
ting energies. The results were significant - even if all 
X-ray emission were synchrotron these Galactic SNRs 
are currently incapable of accelerating electrons beyond 
a limit of 20-100 TeV (for Cassiopeia A the limit is 80 
TeV)^. This work was repeated for 11 SNRs in the Large 
Magellanic Cloud (Hcndrick & Revnolds 2001) this time 
fittin g a Sedov model (Borkowski. Lvcrlv. & Rcvnold^ 
I2OOII) and SRCUT simultaneously, with similar limits 
placed on the ability of SNRs to accelerate electrons. 
These results cast some suspicion on the role of SNRs 
in accelerating ions to energies even below the spectral 
steepening at 10^^ eV in the integrated cosmic-ray spec- 
trum. 

The next step was to move from setting limits to actu- 
ally describing the synchrotron spectrum in SNRs with 
suspected X-ray synchrotron emission. Revnolds (199^ 
showed that the maximum energy attained by electrons 
from the shock acceleration process could be limited by 
several different mechanisms: 1) electrons above some 
energy E^ax could escape from the remnant, e.g. due to 
a lack of MHD waves of the appropriate scale for scat- 
tering, 2) the remnant could be young enough (or small 
enough) that there has not been sufhcient time to ac- 
celerate electrons beyond some Emax, or 3) Emax could 
represent the energy at which radiative losses precisely 
balance gains. The lowest of the energy limits determines 
the nature of the spectrum. These limit s can be ca lcu- 
lated using the expressions given in Revnolds! l|1998() . 

For SN 1006, models limited by rad iative losses wer e 
ruled out even by pre- ASCA data lITievnoldFl IT996|) . 
An estimate of the magnetic field of order 4 /iG 
(from TeV 7-rays as suming an inverse- Compton origin; 
iTanimori et"anll998() . eliminated the age-limited model 
(the magnetic field would have to be below 0.6 fj,G for 
the age limit to be low enough to avoid exceeding the 
observed X-rays). It should be noted that the escape 
model is particularly simple to implement. Unlike the 
loss and age limited models, it is approximately a single- 
parameter model easily adapted for inclusion in X-ray 
spectral analysis software. 

^ IR,evnolds Hr, Keohand I199HI) found that Kes 73 could obtain 
an anomalously high energy of 300 TeV. However, it has since 
been shown that Kes 73 contains a pulsar, i.e. an entirely different 
source of synchrotron emission. 
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In this paper we outline the theory behind the SRESC 
model in fjSland then discuss the SRESC submodels. We 
discuss issues related to X-ray and radio observations in 
21 & ill In we apply the new models to regions of 
SN 1006. In fQwe discuss the results of our fits including 
abundance information, and discuss the implications of 
those results. We then summarize our conclusions in ^Hl 

3. DESCRIPTION OF MODELS 

3.1. Synchrotron model 

The synchrotron models used in this paper consider 
that the shock everywhere accelerates a power-law dis- 
tribution of electrons, with an exponential cutoff above 
an energy i?max whose value depends on various physical 
parameters: remnant age, shock obliquity angle, shock 
speed, and magnetic field strength (details may be found 
in Rcvnolds 19£3). The models then evolve that distri- 
bution behind the shock including adiabatic and radia- 
tive losses, and calculate the volume emissivity of syn- 
chrotron radiation at each point in the remnant by in- 
tegrating that distribution over the single-electron syn- 
chrotron emissivity. Images and total fiux spectra are ob- 
tained by appropriate integrations over the volume emis- 
sivity. All synchrotron models, including SRESC, cap, 
limb, and center, begin with the assumption of Sedov 
dynamics and a power-law distribution of electrons up 
to a maximum energy Emax which in general may vary 
with both physical location within the remnant and with 
time. 

As discussed in we have excluded age and loss 
limited models, in favor of the escape model. In this 
model electrons are presumed to escape upstream, prob- 
ably due to an absence of MHD waves beyond a cer- 
tain wavelength Amax- The wavelength Amax corresponds 
to an energy i?max of electrons. Electrons with gyro- 
radius rg scatter resonantly with waves of wavelengths 
A = 27rrgCOS'0 = 2TT{E/eB) cosip, where ip is the elec- 
tron pitch angle. Therefore electrons will escape up- 
stream once their energy reaches an -Emax given by 

Eniax = Amax 6 Bi/A = 12.0 A17 Bi^G erg, 

where -Bi^ig is the upstream magnetic-field strength in 
units of iiG, A17 is Amax in units of lO^^cm, and we have 
averaged over pitch angles. 

We presume initially that Amax and Bi are uniform 
outside the remnant. Downstream, these electrons radi- 
ate in a magnetic field B2 = rsBi, producing photons 
primarily at a frequency t'max « E^^^B2, more precisely 

z^max - 1.05 X W'^Xl^^Bf^G (^) Hz 

with the magnetic compression ratio varying from 
1, where the shock normal is parallel to the upstream 
magnetic field ("pole"), to the full compression ratio, as- 
sumed to be 4, where the shock is perpendicular (equa- 
torial "belt" ) . So the total spectrum involves a superpo- 
sition of spectra with different turnover frequencies over 
a range of about 4; a fit to an observed spectrum does 
not produce a unique -Emax- This situation is in contrast 
to the much simpler XSPEC model SRCUT, where the 
spectrum is a power-law with an exponential cutoff at 
some Emax- presumed to be the same everywhere in the 
source. 




Fig. 1. — Model and spectrum for (p = 30°. Dashed line is the 
center region, dotted line is the limb regions and dot-dashed line 
is the cap regions. Solid line is the total of all regions. 




Fig. 2. — Model and spectrum for for <p = 60°. Dashed line is 
the center region, dotted line is the limb regions and dot-dashed 
line is the cap regions. Solid line is the total of all regions. 



The implementation in XSPEC of the escape-limited 
synchrotron model SRESC has three parameters: 

1. the radio flux measurement in Janskys at 1 GHz 
(norm) 

2. a, the radio spectral index, where flux density oc 

3. a characteristic roUoff frequency, fitted by XSPEC, 
in Hz (i^roiioff) 

For historical reasons, the frequency froiioff fitted by 
the XSPEC implementation related to fmax above by 

I'rolloff = 5.3 l^max- 

At this frequency the spectrum is about a factor of 6 
below a power law extrapolation from lower frequencies. 
Note that fixing ffoiioff does not determine Xmax, Bi or 
rB independently but only the combination X^^ ^Bf rB- 

The submodel spectra are plotted in Figures [TlSl for 
three values of the aspect angle cj) between Bi and the 
line of sight. The integrated spectrum is insensitive to 
(f), since at these frequencies the emission is extremely 
optically thin, but the submodels differ since the ffux is 
apportioned differently among the regions as a function 
of (j). For instance, the distinction between "cap" and 
"limb" begins to vanish as 4> decreases; at (/) = (not 
shown), the image would be circularly symmetric on the 
sky (though our definitions of cap and limb encompass 
different fractions of the image area). These are spectra 
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Fig. 3. — Model and spectrum for for <j) = 90°. Dashed line is 
the center region, dotted line is the limb regions and dot-dashed 
line is the cap regions. Solid line is the total of all regions. 



for a set of parameters Bi, etc., whose values are arbi- 
trary, for the purpose of comparison among subregions. 
For all values of 0, we notice that even though the limb 
regions cover the smallest area, the limbs still make the 
largest contributions to the total. At low enough frequen- 
cies where the curvature has not yet become important, 
the limb flux is about half the total for all values of (j), 
while the center contribution to the total drops from 0.32 
to 0.19 (as (p drops from 90°) and the cap rises from 0.21 
to 0.34. 

The limb emission is dominated by lines of sight nearly 
tangent to the shell at the very edge, and as the frequency 
increases and the emission is restricted (by electron en- 
ergy losses) to thinner and thinner postshock regions, 
those lines of sight shorten quickly, so that for larger 
values of 4>, the integrated limb spectrum falls off more 
quickly than that of the center, which for larger values of 
(f> also includes the highest-energy electrons (found near 
the "equator"). The softest spectrum is that of the caps, 
which include regions where the obliquity ^Bn is always 
near (parallel shocks); for low values of (j), the center 
does not contain any of the "equator" regions so its spec- 
trum is softer as well. All these differences in the rolloff 
frequencies (the frequency where the decrement from a 
power law is approximately a factor of 6) amount to no 
more than 50% of the input model value of froUoff because 
of the relatively steep dropoff of the spectrum. 

Anticipating the application of these models to 
SN 1006, we note that few objects are likely to be so sym- 
metric that all parts of these models will fit the respective 
regions of the remnants. In the SRCUT and spatially- 
integrated SRESC models already in XSPEC, the models 
are given as decrements below the extrapolation of the 
power-law from radio frequencies. The principal differ- 
ence among the region submodels is in normalization; the 
amount of curvature difference is relatively small, largest 
at ip — 90° and decreasing as (j) decreases. 

In applying these models to SN 1006, we are forced 
to confront a significant problem: as demonstrated in 
Figure 21 at both radio and X-ray wavelengths the model 
overpredicts the emission from the center of SN 1006. 
We will discuss SN 1006's lack of symmetry briefly in 
Section 17^ However, we do not attempt to model the 
detailed radio emission in this paper. We can model 
X-ray emission without presuming a radio structure by 
using the observed radio fluxes and decrements in our 
XSPEC models. That is, we measure the dropoff in the 




Fig. 4. — A comparison of a limb-brightened model (solid line) 
to actual NE-SW slices from the Parkes+Molonglo radio image 
(dash-dot) and the ASCA GIS image (dashes). Note that there is 
not enough central flux in either the X-ray or radio image for the 
morphology of SN 1006 to be explained entirely by limb brighten- 
ing. 



X-ray spectrum from the observed, not the theoretical, 
radio fluxes. 

Similarly, we assume that the same value of Amax char- 
acterizes the emission in all regions. If that is not the 
case, one should in principle produce a totally new the- 
oretical model incorporating some assumed spatial vari- 
ation of Amax- However, given the relatively small vari- 
ations in curvature for different regions, we shall adopt 
a more phenomenological approach, and allow the fit- 
ted rolloff frequencies in the different regions to vary if 
necessary. 

We shall proceed under the assumption that the bi- 
lateral symmetry of SN 1006 is due to an ambient mag- 
netic field close to the plane of the sky, that is, (j) = 90°. 
The differences of model curvatures with obliquity in the 
ASCA band are small enough that this should not be a 
restrictive assumption. In general, substantially longer 
integration times would be necessary for current X-ray 
satellites to detect these subtle differences in curvature. 

3.2. Thermal model 

Since our primary interest is in the nonthermal emis- 
sion, we fit thermal emission from the remnant with 
the simplest reasonable model VPSHOCK. Wherever we 
specify that a thermal model was used, we mean VP- 
SHOCK. This is a plane-parallel shock model with var i- 
able abundances ijBorkowski. Lverlv. & Revnoldsll200lh . 
It represents an improvement over single-temperature, 
single-ionization timescale, non-equilibrium ionization 
models by allowing a distribution of ionization 
timescales. In principle one might expect improvements 
in the results of thermal modeling through the use of 
models at the next level of sophistication, such as Se- 
dov or multitemperature models. However, it is already 
clear from high quality XMM-Newton and Chandra ob- 
servations that thermal spectra demand an even higher 
level of modeling. A promising method is demonstrated 
in Badcnes ct al. ( 2003) where one-dimensional simula- 
tions of different types of explosions are used to gener- 
ate synthetic spectra which are then compared to X-ray 
observations of Tycho. We will limit ourselves to the 
simplest reasonable thermal models, which here play a 
subordinate role to the nonthermal models. 
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4. X-RAY OBSERVATIONS 

Unless stated otherwise SIS & 1 datasets from the 
same observation were fit for each region. Full descrip- 
tions of the datasets used are given in Table ??. We 
used SIS BRIGHT data at high and medium bit rate. 
The standard REV2 screening was used, and data were 
grouped with minimum of 20 counts per channel for valid 
analysis. The background spectra were obtained from 
the 1994 Blank Sky event files in the case of PV 4-chip 
observations and from off source areas for A04 2-chip 
observations. We compared the results of subtracting 
the two types of backgrounds and found the differences 
in the resulting data were minor. What was not minor 
were gain shifts between GIS and SIS data. We intended 
to fit GIS and SIS data simultaneously to obtain better 
statistics. However, we found significant line shifts be- 
tween the GIS and SIS. The difference was 70-80 eV in 
strong line energies (silicon is prominent). Experimen- 
tal GIS matrices reduced the shift in Si line energy by 
~20 eV, but they were unable to completely remove the 
shift. We therefore chose the SIS instrument over the 
GIS for energy reliability and higher spatial and spectral 
resolution. 

Spectra from the GIS do have one advantage over the 
SIS spectra. The GIS detectors have a higher effective 
area at high energy - since this extra signal to noise could 
have implications for the regions dominated by nonther- 
mal emission we verified the results obtained from the 
SIS spectra by checking elliptical GIS regions in the NE 
and southwest (SW). The energy shifts between GIS and 
SIS are less important for measuring the continuum. The 
GIS 2 and GIS 3 data were jointly fit for each limb. We 
obtained radio fluxes (norm) from the same regions used 
to extract X-ray counts, and obtained background from 
regions of similar size outside each respective limb. We 
discuss the results of this comparison at the end of Sec- 
tion 

5. RADIO OBSERVATIONS 

In Paper I, since the entire remnant was being fit, we 
obtained the normalization for the synchrotron models 
from single dish flux measurements. However, in order 
to separate limb from center emission we now needed 
accurate spatially resolved fluxes. Interferometric im- 
ages have the fundamental problem that zero-spacing 
flux (the total flux in the image) is not measured. It 
is well known that maps without this "zero-spacing in- 
formation" are missing the informati on needed to ac - 
curately reconstruct the total flux l)Holdawavl Il999|) . 
The absence of short interferometric spacings showed 
up prominently for SN 1006 as negative flux densi- 
ties in the cent er, NW, and SE regions i n the im- 
ages published inlRevnolds fc Gilmord l|198^ Il993j) and 
iMoffett. Goss fc Reynolds! ()1993|) . ~ero-SDacing infor- 
mation" can only be restored to interferometric maps 
by adding single dish observations. Therefore we ob- 
tained an 843 MHz map created with data from both 
the Molonglo Observatory Synthesis Telescope (MOST), 
an East- West parabolic interferometer with UV cover- 
age from 15 meters to 1.6 kilometers, and the Parkes 
64-meter radio telesco pe, resulting in an image with res- 
olution of 44" X 66" ( Roger et alT ll988l. The total flux 
in the Roger et al. image agrees with single dish mea- 
surements and in this map no regions of the SNR have 



negative fluxes. 

We measured the radio flux in the exact regions from 
which we extracted spectra by using a region file in sky 
coordinates to extract both the spectrum and an X-ray 
image. That image was loaded into AIPS, convolved to 
a resolution of 11" to smooth the X-ray sampling, and 
aligned to the radio image with AIPS task HGEOM. The 
radio image was then clipped everywhere the X-ray image 
was blank, and the remaining fiux was measured at the 
observing frequency (0.834 GHz) and then extrapolated 
to 1 GHz according to the single dish spectral index of 
0.6 (|GreeTill200H) . 

6. FITS TO DATA 

We have fit and plotted the ASCA data from 0.4-10.0 
keV. While 0.4 kcV is below where ASCA data are nor- 
mally considered reliable (generally above 0.6 keV and 
for the most recent data, only above 1.0 keV) we have 
chosen to use these data for several reasons. First, most 
of the datasets used are from the early performance ver- 
ification phase (PV), prior to significant chip damage. 
Second, SISO and SISl are still in good agreement in 
that region - one test for data reliability. Finally, there 
were additional problems fitting regions from the cen- 
ter, southeast (SE) and northwest (NW) in the vicinity 
of 0.6 keV. Without including data from 0.4-0.6 keV, 
XSPEC was free to fit models with extremely high fiux 
at low energies. While we do not trust the data suffi- 
ciently to report abundances measured at low energies, 
we believe the data can be trusted to rule out fiuxes high 
by a factor of two or more. Finally, the lowest energies in 
each fit were dominated by thermal emission. As stated 
above, the main purpose of this work was to test models 
against nonthermal emission, which dominates at higher 
energies. 

There were slight variations reported by SISO and SISl 
in the measured fiux. While could be slightly im- 
proved by allowing the normalizations of the models in 
the two detectors to vary slightly, we did not choose to 
do so, preferring a slightly higher over introducing 
further uncertainties in norm. 

We began by fitting solar-abundance thermal models 
(PSHOCK) to each of the regions (with a limh model in 
the NE and SW). These models fit very poorly, with xl 
ranging from 3.6-20 and obvious line residuals. We can 
say with certainty that even in the bright limbs where 
the nonthermal emission dominates, SN 1006 is not well 
described by solar abundance models. Henceforth, as we 
discuss our fits, all the VPSHOCK models we employ 
allow the individual abundances to vary. 

Table ?? contains the results of the fits to each region 
of the remnant, for each model combination as described 
below. Included are and the number of degrees of 
freedom (DOF, column 1). The following two columns 
determine the nonthermal model: the 1 GHz norm (the 
fiux in Jy, measured from the radio map) in column 3 
(difference between SIS and SIS 1 are a result of gap 
location, rather than chip differences), and in column 4, 
the t'roiioff (^'roiioff in Hz for all SRESC submodels). Ther- 
mal parameters from the VPSHOCK model include: the 
temperature found by each thermal model in units of keV 
(column 5) and the ionization timescale r = riet in units 
of s cm~^ (column 6). In column 7 we give the emis- 
sion measure (norm), defined as ^fjyr J neUndV where 



6 



Dyer et al. 




Fig. 5.— Thermal and nonthermal fits to the northeast limb. FiG. 7.— The northeast and southwest limbs with the maximum 

amount of nonthermal synchrotron plotted. 




Fig. 6. — Thermal and nonthermal fits to the southwest limb. 



D is the distance to the source in centimeters. The fitted 
abundances are listed in column 8-12, given as ^^^"^^ ^ 

(by number, not mass, from iGrevesse fc Andera 11983) . 
Since the hydrogen abundance is poorly constrained by 
the fits, abundances were fit relative to silicon. Errors 
given in Table ?? are l.Sa errors (87% confidence in- 
terval). As in Paper I, for both GIS and SIS data, we 
assumed an absorption column density of 5x10^° cm~^, 
using the Wisconsin absorption model, wahs. 

Note that, except where specified, all synchrotron mod- 
els were fit with normalization and spectral index fixed, 
as measured from radio observations. There is consid- 
erable degeneracy in the model between norm, a, and 
z/roiioff: demonstrated in Table ??, but radio measure- 
ments allow us to input, rather than fit, norm & a. 

6.1. Brighter limbs: northeast and southwest 

We began by searching for firm upper bounds to the 
synchrotron emission from the bright limb regions, shown 
in Figures |S1 & We did this by finding the maximum 
amount of the limb model tolerated by the SIS NE and 
SW regions, assuming a spectral index of 0.60 and limb 
norm measured from the radio, i.e. we raised t'roUoff 
until it threatened to exceed the data at high energies. 
We found that the spectra of the NE and SW regions 
are not identical (see FigureEj). The NE could tolerate a 
maximum t'roUoff of 7x10^^ Hz while the SW could only 
tolerate a maximum t'roUoff of 4x10^^" Hz. 

We then fit the data with combined thermal and 
nonthermal models, using VPSHOCK (variable abun- 
dances) -I- /im6, shown in Figure |S| The results for the 



fit are given in Table ??, row 1. The combined thermal 
and nonthermal fit achieved a xj~1.8. Residuals from 
the fits primarily show discrepancies between SISO and 
SISl between 0.5 and 1.0 keV. 

With the addition of a thermal model the froiioff for the 
NE dropped to 3.3x lO^"^ Hz with 1.5a limits of 2.6x lO^"^ 
and 3.8x10" Hz. The SW had a u.oiioS of 2.3x10^^ Hz 
with 1.5(7 limits of 1.9x10" and 2.6x10" Hz. These 
fitted values of t'roUoff do appear significantly different, 
though not by as much as the maximum values above. 
This agrees with the higher radio-to-X-ray ratio in the 
SE shown in Figure 0] In order for the radio-to-X-ray 
ratio to be lower than in the NE, ^'roiioff' must be lower 
(since a is fixed.) The thermal component of the NE had 
a temperature of 1.85} gf keV and ionization timescale 
of 5.2|:| X 10^ s cm-3.' The total flux from 0.4 to 10.0 
keV in the NE was 6.1x10^^^ ergs cm~^ s^^, of which 
47% was from the synchrotron model. The SW had a 
temperature of 1.87j 7g keV and an ionization timescale 
of 1.5};! X lO^^cm^^' s. The total flux in the SW was 
4.3x10"^^ ergs cm~^ s~^, of which 46% was from the 
synchrotron model. 

The GIS data, in contrast to the SIS data, show no ob- 
vious lines between 2 and 10 keV. Therefore there is lit- 
tle purchase for a thermal -f nonthermal fit, although for 
comparison with SIS results we did try fits with SRESC 
and SRESC+VPSHOCK models. Since the GIS has a 
much wider bandpass than the SIS we also allowed the 
spectral index to vary, in order to eliminate the possi- 
bility that freezing a within the narrower SIS bandpass 
could introduce false results in froiioff (as discussed above 
and shown in Table ?? the degeneracy in norm, ^'roiioff, 
and a decreases the usefulness of fitting a). 

The 1 GHz norms from the GIS regions were 3.37 Jy 
for the NE and 5.11 Jy for the SW. (AU errors on GIS 
observations are 1.65 a, i.e. 90% confidence interval.) 
With a fixed to 0.6, we obtained values of I'roiioff of 
5.89|:f| X 10" Hz in the NE and I'l-oUoff = 3.25i:?[] x 10" 
Hz in the SW. With a allowed to vary we obtained 
= 2.6li ?| X 10" Hz and a = 0.55 ± 0.01 in the 
NE, and z/roiioff = 2.62^:^i x lO" Hz and a = 0.59±0.01 in 
the SW. The thermal-|-nonthermal fits were carried out 
with SRESC plus a solar abundance VPSHOCK model. 
With a fixed at 0.6, we found t^roiioff to be 4.503 57 ^ ^^^^ 
Hz in the NE and 2.65?:S]| x lO" Hz in the SW. 

This confirms results found from the SIS fits, includ- 
ing the slightly higher roUoff frequency in the NE. In the 
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Fig. 8. — Northwest: a) SIS region shown on the Parkes radio 
image b) high cap model with i^rollofl measured at the limbs c) 
thermal model VPSHOCK. 





Fig. 10. — Center: a) SIS region shown on the Parkes radio 
image b) center model emission predicted by the i^rollofl measured 
at the Umbs c) thermal model VPSHOCK. 




Fig. 9. — Southeast: a) SIS region shown on the Parkes radio 
image b) cap model emission predicted by the t'rollofi measured at 
the Umbs c) thermal model VPSHOCK. 



thermal+nonthermal fits the GIS do not show any signif- 
icant differences in t'roUoff that might be expected from 
the better GIS sensitivity at high energies. The results 
are consistent with the results for SIS fits, which are pre- 
ferred since they more accurately constrain the thermal 
emission. 

6.2. Fainter limbs: northwest and southeast 
Northwest 

Our initial assumption in applying the simple escape 
model was that all regions of the remnant would be fit 
with varying amounts of norm, measured in radio ob- 
servations, but the same spectral index and fc'roiioff- We 
note that this procedure normalizes away any departure 
of the radio morphology from model predictions (for ex- 
ample, the simple SRESC model predicts higher bright- 
ness in the remnant center than is observed). It sim- 
ply requires that no additional departures from model 
predictions would be necessary between radio and X-ray 
wavelengths. 

However, it is immediately clear that while the NE 
and SW regions have similar t'roiioff's (average 2.8x10^'' 
Hz), that this value of i/roUoff in the cap submodel, with 
norm appropriate to the NW and SE, clearly contradicts 
the NW and SE data, even before addition of a ther- 
mal model required by the presence of lines (see Figures 
IHh) and Eb) . This agrees wi th the high spatial resolu- 
tion XMM-Newton results of iDecourchelld l|2002D which 
found ^^roiioff varying with both radius and azimuth. 

Fits to the NW region with a thermal-only variable 
abundance model (shown in Figure O;) yielded a xt of 
1.80. Deficiencies in the thermal model can be noted in 
residuals near 0.7 and 0.9 keV (Fe L-shell or Ne Ka). Full 
results are given in Table ??, row 5. The temperature 
given by the fit was 1.2} q keV and r — G.Sg ? x 10^ cm~'^ 
s. The total flux from this fit was 1.3xl0~ ergs cm~^ 
s""'^. While a fit with thermal+nonthermal models was 
attempted the data preferred no nonthermal component. 

It was critical to place upper limits on the amount 



of nonthermal emission even in regions where no non- 
thermal component was required so we used the follow- 
ing procedure: we added a nonthermal model (beginning 
with a = 0.60, ^'roiioff=2.8x 10^^ Hz and norms measured 
from the radio image) to the best thermal-only fit (given 
in Table ??, row 5) and then changed each of the three 
nonthermal parameters in turn (a, i^roiioffj and norm), 
holding the others fixed until the rose by ^--^2.3, a 
1.5a change. The results of these fits are given in Ta- 
ble ?? with the varying quantity in bold face. With this 
method we found that for the NW data the cap model 
had a maximum t'roUoff of l.lxlO^'' Hz, or a maximum 
norm of 0.25 Jy at 1 GHz, or a spectral index a of 0.67 
or steeper. This places a limit of <8% on the nonthermal 
contribution to the NW flux. 

Southeast 

As with the NW, using the value of ^'loiioff measured 
in the NE and SW regions overpredicted the observed 
flux in the SE (see Figure Eb). Thermal fits with vari- 
able abundances to the SE (see Figure IHli) were signifi- 
cantly worse than the NW. xt was ~3 for a single VP- 
SHOCK model, a two- VPSHOCK model (not shown), 
and a capH- VPSHOCK model. Interestingly, despite a 
much worse xt the residuals look very similar to those in 
the NE - obvious line-like residuals at 0.7 and 0.9 keV. 

The thermal-only model (Table ??, row 7) had a tem- 
perature of l.lJ o keV and a t — 2.3f | x 10^" cm^^ s. 
As above for the NW, while a combined thermal and 
nonthermal model was tried, the data preferred no non- 
thermal component at all. 

Following the method outlined above for the NW, we 
set limits on the nonthermal emission. In the SE we can 
exclude a ^'roiioff above l.lxlO^'' Hz or a spectral index, 
a fiatter than 0.67, or a norm above 0.2 Jy at 1 GHz. 
The total fiux of this faint region is 8.6x 10^^^ ergs cm~^ 
s^^ and the maximum nonthermal contribution is < 9%. 

6.3. Center 

We turn now to the center of SN 1006. Unlike the 
NW and SE, the center data could just tolerate the non- 
thermal model center with fixed t'loUoff , shown in Figure 
llOb . However, joint fits soon reduced the nonthermal 
component to <5% (using the method outlined in fits 
for the NW and SE). Both combined and thermal-only 
models had very high xt{> 5). Thermal-only fits (Fig- 
ure and Table ??, row 3) reported a temperature of 
0.82[j:?| keV and a r = 2.5^:| x 10^° cni^^ s. As with the 
NW and SE, when a thermal+nonthermal model fit was 
attempted the data preferred no nonthermal component. 
All fits to the center show residuals similar to the NW 
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and SE. 

The thermal fits for each region of the SNR reveal ele- 
vated abundances. We can make a separate case for su- 
persolar abundances by measuring the equivalent width 
of the well-separated Si Ka line in the center where the 
lines are strongest. The result is an equivalent width for 
silicon of 600 ± 50 eV (1.65 ct, 90% confidence). We draw 
conclusions from this in Scction l7.2l 

7. DISCUSSION 
7.1. Difference between the NE and SW limbs 

The NE and SW regions do have slightly different spec- 
tra, with values of i^ioUoff differing by a significant amount 
in both SIS and GIS fits including a thermal component. 
(The froUoff values bracket the Paper I value averaged 
over the remnant of 3.1 x 10^^ Hz.) This is not, however, 
sufficient to explain the differences in TeV detection of 
the NE and SW regions (assuming as in Paper I that 
the TeV detection is cosmic microwave background pho- 
tons, inverse-Compton scattered by the same population 
of electrons that produce the X-ray synchrotron emis- 
sion). The small difference in t'roUoff between NE and 
SW regions only produces a difference of 5% in this flux 
at 1 TeV. Even t'roUoff's differing by factor of three only 
produce a difference of VS ^ 70% in E„iax- A differ- 
ence in the spectral index is far more likely to produce 
the differences between synchrotron emission in the NE 
and SW regions which would be required to explain dif- 
ferences in 7-ray detection (although, as discussed below 
the spectral index cannot account for the entire change). 

Using the formula given in t l3.1l and assuming a mag- 
netic compression ratio r=4 and value Bi— 3 fiG as 
found in Paper I, we find that the I'roiioff of S.SfJ x 10^'' 
Hz in the NE gives a Amaa;=l-5i'3 x 10"'^^ cm and in the 
SW a froiioff of I X 10" Hz gives Xmax^l-^l} x 10^^ 
cm. These values imply electron escape above 33 TeV in 
the NE or 28 TeV in the SW, similar to the value of 32 
TeV quoted in Paper I, as is to be expected. 

7.2. Thermal fits 

As stated earlier, the focus of this work was on testing 
nonthermal models. Therefore only the simplest plane 
shock model was used to fit thermal emission. We ruled 
out solar-abundance thermal emission in each region of 
the SNR and then used the variable-abundance model, 
VPSHOCK, for the remainder of our fits. 

The quality of the fits varied widely across the SNR. 
In bright regions such as the NE, the SW, and the NW, 
the fits gave marginal values of xt ~ 1-8. Fainter regions 
such as the center and the SE were poorly fit by VP- 
SHOCK with xi > 5.0. However, despite different values 
of xtj the residuals of all regions with significant thermal 
emission (SE, NW and center) had line-like residuals at 
0.7 and 0.9 keV. This casts suspicion on even the better 
fits and suggests problems with the underlying atomic 
data since this is where 3s to 2p transitions in Ne-like Fe 
ions are found. These transitions appear to be stronger 
than expe cted in astrophysical sources , including SNRs 
(as seen in Ivan der Hevden et al.ll2'002() . 

In the current NEI vl.l models in XSPEC which 
we use, Fe L-shell data are based o n calculations by 
iLiedahl. Osterheld. fc GoldsteiiJ l)1995i|l for electron col- 
lisional excitation of Ne- to Li-like Fe ions, and on older 



atomic data for inner-shell processes in Mg-like and Na- 
like Fe ions l|Hamilton fc Sarazinlll98^ . Indirect pro- 
cesses such as inner-shell collisional ionization, radiative 
and dielectronic recombination, and resonance excitation 
are not included, although they have been shown to be 
important in enhanci ng 3s to 2p transitions in the Ne- 
like ion in particular l)Gull200^ . The enhancement of 3s 
to 2p transitions in the Ne-like ion due to the inner-shell 
ionization of the Na-like ion might be particularly im- 
portant in high-temperature plasmas at low ionization 
ages, such as those en countered in SN 1006. In addi- 
tion, EEi^OD l|2003t) showed recently that atomic res- 
onances in the electron collisional excitation of Ne-like Fe 
are also important, although they have been neglected 
in the past. So there are many deficiencies in Fe L-shell 
atomic data which could result in poor fits in the 0.7-0.9 
keV range. 

We tested this hypothesis by adding lines arising from 
the inner-shell ionization of the Na-like Fe ion into 
XSPEC NEI v2.0 models. These models are based on 
the Astrophysica l Plasma Emission Database (APED; 
iSmith et alJ2001|) . and are currently missing several pro- 
cesses important in NEI conditions, but they allow for 
easy modification of the atomic database. We added 
three lines of the Ne-like ion at 17.098, 17.053, and 16.778 
A (0.725, 0.727, and 0.739 keV, the strongest at 0.725 
keV). For a shock model with parameters appropriate 
for the NW region of SN 1006 (fcT = 1.17 keV and 
T = 6.5 X lO^cm"'^), their contribution to the Fe spec- 
trum is substantial. So it appears certain that the cur- 
rent fits to SN 1006 spectra use inadequate Fe L-shcU 
atomic data, and that one problem with current ther- 
mal fits can be traced to Fe L-shell data. In particu- 
lar, missing lines from inner-shell ionization of Na-like 
Fe ions are certainly a problem for the NW, and most 
likely elsewhere. Since this paper is primarily concerned 
with the nonthermal emission and better spectral infor- 
mation (for smaller spatial scales) has been obtained by 
XMM-Newton and Chandra we defer advanced models 
to a future paper. 

This may not be the only issue preventing us from 
finding adequate the rmal fits. X-ray sp ectral model- 
ing of Type la SNRs l|Badenes et alJl2003D has revealed 
significant differences in physical conditions between 
chemically-distinct layers of shocked ejecta. In partic- 
ular, the emission-measure averaged ionization timescale 
of Fe appears to be shorter than that of other abundant 
elemen ts in most models considered by iBadenes et alJ 
l)2003l Fig. 6). It is then unlikely that the chemically 
homogeneous VPSHOCK model can provide a satisfac- 
tory fit to both Fe L-shell lines and lines produced by 
other elements. Substantial improvements in both Fe L- 
shell atomic data and in spectral models are required in 
order to make further progress in modeling SN ejecta in 
SN 1006. 

The Paper I fits to the total flux of SN 1006 achieved 
better values, x^=l-2. However, the better fits are causes 
by two factors: the poorer spectral resolution of the 
ASCA GIS instrument and the fact that the integrated 
spectrum is dominated by the synchrotron contributions 
of the two bright limbs which are well described by our 
models. 

One conclusion that can be drawn from the thermal 
fits in Table ?? is that while the thermal fits may be 
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oversimplified, the VPSHOCK model is clearly sensitive 
to both thermal parameters and abundances. Whereas 
we might have expected large uncertainties, the 1.5a sta- 
tistical errors are rarely more than 8%. Clearly, the data 
are capable of constraining thermal models; the correct 
thermal model should yield significant information. 

While our thermal model fits to the center were of poor 
quality statistically, they strongly implied elevated abun- 
dances of heavy elements such as Si. Fortunately, a more 
model-independent criterion is available: the equivalent 
width of the Si Ka complex (Si XIII). Our measured 
value of 600 ± 50 eV is higher than a plane shock with 
solar abundances could achieve until r ~ 3 x 10^'^ cm~'^ 
s, by which time the shock has kT considerably above 
1 keV, and should have much hi gher ratios of O VII to 
Si XIII than we observe ijLong et a l. 2003). We conclude 
that the central emission is dominated by shocked ejecta. 

7.3. Nonthermal results 

Our spatially resolved subsets of SRESC, cap, limb, 
and center, differ surprisingly little from one another or 
from the total (SRESC itself), though the differences we 
do see (' tl3.1|l are in the sense one would expect. How- 
ever, for these ASCA data, and in fact for most data from 
current X-ray telescopes, the differences between the spa- 
tially integrated SRESC model and the submodels are 
not likely to be distinguishable. This could not have been 
predicted in advance, but it is fortunate for data analysis. 
However, our conclusions about the lack of cylindrical 
symmetry are stronger for our having based them on the 
appropriate submodels that symmetry would require. 

The most striking result from the spatially resolved 
fit is that the SRESC submodels with values of froUoS 
that fit the limbs, do not fit the spectrum in the cen- 
ter, NW and SE. There are several possible explanations, 
with interesting consequences. First, it is possible that 
the SRESC submodels are an accurate description of the 
nonthermal emission and that the value of J^roUoff differs 
from one region to the next. Since I'roiioff oc X^B^ a 
difference in fc'roiioff would imply a change in the exter- 
nal magnetic field or the maximum wavelength of MHD 
waves available for scattering. The variation of the t'roUoff 
parameter across the remnant marks a significant depar- 
ture from the simple assumptions of the escape model. 
We can offer no justification for this but if B or Xmax 
could be isolated they would be important probes for 
the shock physics. 

As demonstrated in Table ??, the SRESC submod- 
els show some degeneracy between J^roUoff and the radio 
spectral index a, in that a given X-ray to radio flux ratio 
can be produced (within limits) by flatter spectra with 
lower froiioff, or steeper spectra with higher froiioff- A 
large source of error in all synchrotron fits is the uncer- 
tainty in a. An uncertainty in a of ± 0.1 at 1 GHz can 
change the extrapolated flux in the 2 to 3 keV range by 
a factor of 7. In addition to being poorly measured by 
single dish radio observations, the spectral index could 
vary across the source. We can, however, rule out the 
possibility that spectral index is solely responsible for 
the variation in froiioff- While radio observations cannot 
constrain a to better than ±0.02, single-dish observa- 
tions flrmly exclude an a as steep as 0.67. Since the two 
limbs dominate the total radio emission as they do the 
X-ray emission, the kind of difference between limbs that 



would be required seems unlikely. We cannot, of course, 
rule out a contrived situation in which the spectrum of 
the SW suddenly steepens significantly compared to the 
NE at radio frequencies higher than have been measured. 
Nevertheless, the question of a varying spectral index 
across the remnant is intriguing. Many researchers have 
searched for variations in radio spectral index in SNRs 
as an indication of changes in shock acceleration. How- 
ever, interferometric observatio ns are ill suited to this 
kind of study (see discussion in Dver fc R evnolds 1999) 
and current observations of SN 1006 can neither support 
nor contradict this premise. 

A second possibility is that the physics in the syn- 
chrotron model that underlies SRESC and submodels is 
inadequate to describe the nonthermal emission. One 
critical assumption is that the diffusion coefficient n is 
proportional to the energy. A divergence from this as- 
sumption at the highest energies could have serious con- 
sequences. 

7.4. Morphology 

We began our investigation by assuming cylindrical 
symmetry about the NW-SE axis. We can now rule this 
out with a fair amount of certainty. If SN 1006 were 
cylindrically symmetric a line of sight through the center 
would sample a predictable amount of the nonthermal 
emission which brightens the limbs. This is ruled out by 
the fact that the SRESC submodels over-predict the high 
energy flux in the NW, SE and center, as demonstrated 
in Figures|SjD,|Hb andllOb. This deficit of central emission 
is in addition to the deficit already present between ob- 
served and theoretical radio images in the SRESC model, 
since we are using measured rather than theoretical val- 
ues of norm. An explanation which might account for the 
weaker central emission in radio compared to the simple 
model (for instance a substantial radial component to the 
magnetic field) cannot explain the additional deficit in X- 
rays. Any explanation must impact the X-ray emission 
differently than the radio emission. 

That SN 1006 is not simply limb brightened due to a 
NW-SE symmetry is an uncomfortable conclusion, since 
it implies a preferred orientation between SN 1006 and 
the observer, al though there ar e othe r arguments for 
asymmetry from lWillingale et al.l l)1996() . based on ratios 
between the X-ray emission in limbs and center as ob- 
served by ROSAT, and arguments derive d from the spec- 
tra o f the Schweizer-Middleditch star ijHamilton et alJ 

This asymmetry could take several forms. It is pos- 
sible, likely in fact, that the upstream magnetic field is 
not uniform. There is already evidence for asymmetry 
in the NW, where Ha is strongest - presumably the am- 
bient medium outside the NW is, for some reason, less 
ionized. This would have an impact on electron acceler- 
ation. This could affect Xmax in particular. 

It is possible, of course, that fundamental assumptions 
made in the SRESC model are incorrect. While the 
agreement of the model with the bright limbs argues that 
the basic picture of electron acceleration to a power-law 
with an exponential cutoff is sound, at least in those re- 
gions, it is possible that a completely different geometry 
holds, for instance, one in which the bright limbs are 
"caps" seen edge-on rather than brightened limbs of an 
equatorial "belt." The symmetry axis would then run 
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NE-SW, but would need to be nearly in the plane of 
the sky. The reasons that the polar caps are brighter 
in synchrotron radiation in such a picture would need to 
involve some superior feature of parallel shocks for elec- 
tron injection or acceleration, or perhaps magnetic -field 
amplification as sugg ested bv iBell fc Lilce^ l^MH^ . A 
fundamental defect of such a picture is that remnants for 
which the ambient magnetic field is closer to the line of 
sight should appear in radio synchrotron images as two 
maxima surrounded by a steep-spectrum halo, a mor- 
phology which is not seen in any Galactic remnant. 

7.5. Comparison with Chandra Results 

Chandra observed part of the NE limb in Cycle 1 and 
part of the NW limb in Cycle 2. The smaller field of view, 
and charge-transfer-inefficiency problems for the front- 
illuminated CCDs, restricted initial analysis to much 
smaller regi ons than those w e analyzed here. However, 
the results ijLong et al.l 12003(1 are consistent with ours. 
In the NE, thin filaments appeared to have completely 
line-free spectra, and were well described by an SRESC 
model with a = 0.54 and froUoff of 6.9 x 10^^ Hz. Since 
Chandra was able to isolate a region immediately be- 
hind the shock only a few tens of arcsec in size, it is to 
be expected that the hardest electron spectrum (high- 
est froiioff) would be found there. Our CIS value for 
froiioff, with a allowed to float, is close to this value, al- 
though slightly lower because of dilution from thermal 
emission and nonthermal emission further downstream 
where adiabatic and radiative losses begin to soften the 
spectrum. A region several arcminutes behind the shock 
in the NE was shown to have substantial line emission, 
and was fairly well described by a VPSHOCK+SRESC 
model with fcT - 1.4 keV, r - 3 x 10^ cm^^ j^j^j 
a t'roiioff of (3 — 6) X 10^^ Hz. Elevated abundances of 
heavy elements (chiefly Si and Fe) were required. These 
results are similar to our inferences about the thermal 
component in the NE, though these of course apply to 
a much larger region than the Chandra results. In the 
NW, a VPSHOCK fit to a small region coincident with 
the Ha optical filament gave kT ~ 0.9 keV, r ~ 6 x 10^ 
cm^^ s, and required no nonthermal component (though 
it could tolerate a weak one). These results are compa- 
rable to those we report here from a much larger region. 
The broad similarities between fitted values of small re- 
gions with Chandra data and larger averages with ASCA 
data give some confidence that the systematic residu- 
als we observe that degrade the statistical quality of our 
fits are not spacecraft-dependent systematic effects. (Of 
course, they may be inherent in the models or atomic 
data, which were the same in the two analyses.) While 
Chandra observations did include regions closer to the 
center of SN 1006, those were on ACIS-I chips subject to 
charge-transfer-inefRciency problems, and the data have 
not yet been analyzed. 

7.6. Effect on Paper I Results 

These new investigations change few of the results in 
Paper I. It is clear that synchrotron emission from the 
limbs dominates the total fit: the value of froiioff in Paper 
I was 3.02 8 ^ 10^^ Hz, while fits to the limbs individually 
found 3.3 X 10^^ and 2.3 x 10^^ Hz. 

While in theory more precise values for the postshock 
magnetic field strength, the energy in relativistic elec- 



trons and the electron efhciency (9 /^G, TxlO**^ ergs and 
5% in Paper I) could be obtained from limb-only fits, 
in fact these will change insignificantly considering other 
uncertainties. 

8. CONCLUSIONS 
8.1. Nonthermal Emission 

1. From comparing the limb-to-center ratio in X-ray 
and radio wavelengths (Figure^, we find that the mor- 
phology of SN 1006 is not solely due to limb brightening 
of a cylindrically symmetric source about a line parallel 
to the bright limbs. This has the inescapable consequence 
that the observer is now required to have a particular 
vantage point with respect to the SNR. 

2. SRESC provides a useful framework to describe 
synchrotron emission in SNR blast waves. Spatially re- 
solved subsets of the SRESC model, cap, limb, and cen- 
ter, differ surprisingly little from the spatially integrated 
model, with the greatest differences occurring when the 
upstream magnetic field is in the plane of the sky. The 
differences are unlikely to be detectable with the current 
generation of X-ray telescopes. 

3. The limb model provides a good description of the 
limbs of SN 1006, but the cap and center region models 
overpredict nonthermal emission in other regions of the 
remnant if I'roiioff is assumed to be constant everywhere. 

4. We do find differences in the nonthermal spectrum 
between the NE and SW, but they are not sufficient to 
explain why only the NE is detected in TeV 7-rays, as- 
suming the 7-rays are produced by inverse-Compton up- 
scattering of the cosmic microwave background photons. 

The usefulness of limb, cap, and center models for 
other SNRs is not obvious. Most SNRs have more 
complex morphology than SN 1006, and few present 
regions clearly appropriate for a center or cap model. 
While the limb model may be appropriate for filamen- 
tary structure at the edge of other SNRs, its curva- 
ture differs only very slightly from the SRESC model. 
We recall that the model SRESC, and its derivatives 
have several assumptions that make them unsuitable 
for SNRs of unknown type - the models presume a 
Sedov-phase remnant encountering a uniform upstream 
medium with a constant magnet ic field. SRCUT, as dis- 
cussed in Revnolds fc Keohand l |l999.) , which describes 
synchrotron radiation with a minimum of assumptions, 
may be a better choice for SNRs about whose environ- 
ment less is known. 

8.2. Thermal Emission 

1. We can rule out solar abundances in the thermal 
models in all regions of the SNR, including regions dom- 
inated by nonthermal emission. 

2. Simple thermal models are incapable of describ- 
ing the ASCA data adequately. We note that in much 
smaller regions of SN 1006, th e analysis of Chan dra data 
encountered similar problems l|Long et alJl2003i) . We be- 
lieve that this is mainly caused by deficient Fe L-shell 
atomic data, where many important processes (such as 
inner-shell ionization of the Na-like Fe ion) are currently 
not included in NEI spectral codes. It is also likely that 
the available spectral models are not adequate for model- 
ing chemically-inhomogeneous SN ejecta. A full analysis 
of the combined thermal and nonthermal emission from 
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SN 1006 will require XMM-Newton and Chandra quality 
data with substantial improvements in atomic data and 
modeling. 

3. To the extent that our thermal results can be 
trusted, we see elevated abundances which suggest the 
thermal emission contains substantial contributions from 
the reverse shock. Here the models agree with the sili- 
con equivalent width, and with the results of other X-ray 
analyses of SN 1006. 
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Table 1: Datasets Fit For This Paper 



Observation Date Cycle t [ks] ceiiter[a, S] region size ['] Count rate; SISO SISl 

NB Limb 1993Septl3 PV 21 15 03 45, -41 47 00 23x9 2.38/2.33 2.03/1.96 

SWLimb 1993Augl9 PV 26 15 02 00,-42 00 15x8 1.47/1.42 1.23/1.18 

Center 1993Augl9 PV 26 15 03 00,-41 56 00 16x13 1.04/0.98 0.91/0.83 

NW 1993Augl9 PV 26 15 02 15, -41 35 46 8x8 0.25/0.22 0.30/0.27 

SB 1996Feb21 A04 26 15 03 15, -42 05 30 4x8 0.19/0.18 0.18/0.16 
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Table 2. Best Fit Parameters to ASCA SIS Observations of SN 1006 



Nonthermal SRESC Thermal VPSHOCK 



Region 



kT [keV] 



NE Limb |S| 4.60,4.84 3.3E17|;||;J5 



1.85j 
1.87} 



5.2E9|:|g| 1.5E-2;;J2_^ l.OE-1^ 



4.1E-1° 
1.7E-lJ 



4.6E-1^ 
5.8E-l5 



I.IE-2I 



0.82q 11 2.5E102;|eio 2.9E-22;7g_2 l-12E-lj' jgg_ J 5.5E-2g;jg_2 1.8E-lj °g_j 8.6E-lg ^ 4.4E-2^ ! 



l.lEOl;||"_j 5.9E-2^;; 



"^Measured fluxes are given for SISO and SISl. The difference is primarily due to the location of chip gaps. 
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Table 3. Upper Limits on SRESC Models in SN 1006 



Region 


Model 


a 


l'rolloff[Hz] 


norm, [Jy@l GHz] 


Northwest 


cap 


0.60 


2.8x10" 


0.247 






0.60 


1.068x10" 


1.01/0.93 






0.669 


2.8x10" 


1.01/0.93 


Southeast 


cap 


0.60 


2.8x10" 


0.192 






0.60 


1.023x10" 


0.76/0.95 






0.671 


2.8x10" 


0.76/0.95 


Center 


center 


0.60 


2.8x10" 


0.42 






0.60 


6.9x10''^ 


2.51/2.57 






0.692 


2.8x10" 


2.51/2.57 



Table 4. Percentage Thermal and Nonthermal Emission in Regions of SN 1006 



Region 


Thermal flux 


NT Flux 


Percent NT 




10"" ergs cm"^ s"^ 


10"^^ ergs cm"^ s"^ 




NE Limb 


3.24 


2.83 


47% 


SW Limb 


2.33 


1.96 


46% 


Center 


3.27 


<0.18 


<5% 


NW 


1.25 


<0.11 


<8% 


SB 


0.86 


<0.09 


<9% 



Fig. 13.— 



